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Abstract 



The dynamics of the pion field after a quench is studied in the frame- 

£> | work of the linear sigma model. Our aim is to determine to what extent 

E***- ' the amplified pion field resembles the DCC picture originally proposed 

,__i . in the early '90s. We present the result of a computer experiment where, 

0^ \ among other things, we study in detail the correlation between isospin 

orientations of the distinct modes of the field. We show that this corre- 
(^ ■ lation is absent. In a sense, the distinct modes behave as distinct DCCs. 

"*^h ! The implications of this observation are discussed. 

6-t, 

1 Introduction 

X 

A disoriented chiral condensate (DCC) is a medium where "the quark con- 
densate, (0|gLgij|0), is chirally rotated from its usual orientation in isospin 
space" 0]. This hypothetical object has attracted the attention of many physi- 
cists in the last decade and is still a topic of intense theoretical as well as 
experimental investigations (see reviews and also ||). 

One imagines a rapidly expanding fireball, produced in a high energy 
hadronic or nuclear collision, whose interior is separated from the outside vac- 
uum by a thin, hot shell of hadronic debris. The quark condensate inside 
the bubble might then be misaligned with respect to its vacuum orientation. 
The DCC subsequently decays toward ordinary vacuum by radiation of soft 
coherent pions. The original picture of the DCC is that of a classical pion field 
oscillating coherently in a well defined direction in isospin space [|J and indeed, 
Blaizot and Krzywicki showed that such field configurations can develop with 
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suitable initial conditions ||. This phenomenological picture leads in partic- 
ular to a very simple, but very striking prediction: the 1/a/7 event-by-event 
distribution of the neutral fraction / of radiated pions. 

The idea received a boost in 1993 when Rajagopal and Wilczek pro- 
posed a microscopic mechanism for the formation of a classical coherent pion 
field inside the bubble 0, the so-called quench scenario: the rapid expansion 
causes the system, initially thermalised above the critical temperature, to cool 
down very rapidly, leading to the amplification of the soft pion modes via 
the mechanism of spinodal instability. Since then, this scenario has attracted 
much attention and has been further developed. In particular, quantum effects 
in the mean field approximation have been included (see e.g. 0, ||) and the 
drastic quench approximation has been abandoned in favour of one taking the 
expansion explicitly into account || U . 

The quench scenario has been widely accepted as a microscopic descrip- 
tion of DCC formation in heavy ion collisions and some works were dedicated 
to the study of its phenomenological implications [|TI| [TT|, |T2|. However it is 
less clear than usually believed whether the typical field configuration emerg- 
ing from a quenched thermal ensemble is identical to what one had in mind 
in the early '90s. Indeed the DCC configuration can be characterised by three 
essential features QTJ, §, f|: 

(a) it is a coherent field excitation, essentially a classical field configuration, 

(b) each Fourrier component of the field oscillates along a well defined orien- 
tation in isospin space?], 

(c) the orientations of different modes are aligned along one single direction. 

The object of the present paper is to perform a detailed statistical analysis of 
the generic field configuration emerging from the quench scenario in order to 
determine whether it exhibits a DCC structure. 

For this purpose one has to choose a reliable framework. Quantum effects 
are not of first importance as the system is essentially classical || . Moreover 
these corrections can, so far, only be included in a mean field approximation, 
not suitable to study the correlations between modes. Finally, in a realistic 
scenario, where the quench is due to a rapid expansion, the system enters only 
rarely in the instability region [[13]. Then, to keep the argument simple, it 
is worth going back to the original framework of Rajagopal and Wilczek in 
which the system, initially in a high temperature (T > T c ) configuration is 
evolved with the full, non-linear, zero temperature equations of the classical 
linear a-model. 



2 In the following we refer to the trajectory of one Fourier mode in isospin space as 
(iso-)polarisation. The case described in (b) will then be refered to as linear polarisation. 



2 The formalism 

2.1 The quench scenario 

The chiral field is parametrized by a four-component vector field in chiral space: 
4> = (-7T, er) and its time evolution is governed by the classical linear a-model 
equations of motion 

(d 2 - Xv 2 + \(f> 2 (x, £)) <j>(x , t) = Hu a , (1) 

where the parameters v, A and H are related to physical quantities via: 

2Xf 2 = X (f 2 - v 2 ) , and H = U 



ml = ml - 2Xfl = A {ft - v 2 ) , and tl = Um~^ , 



and are chosen so that m 7r =135 MeV, m CT =600 MeV and / 7r =92.5 MeV. 

The initial field configuration is sampled from a thermal ensemble at 
some temperature T > T c : the values of <f> and (ft are chosen independently 
on each site of the cubic lattice from gaussian distributions centered around 
= = and with variances^ (<f) 2 ) = v 2 : /16 and (^) = v 2 /A (j = 1, 2, 3, 4). 
The lattice spacing a has then the physical interpretation of the degenerate 7r 
and a correlation length at temperature T, a = (200 MeV) -1 . 

We denote by <p(k,t) and ip(k,t) the Fourier components of the field 
and its time derivative at time t, and choose Neumann boundary conditions 
so that these components are real. This is not an essential point, but it is 
more convenient to make this choice in discussing the question of polarization. 
The infrared cutoff in Fourier space is Afc = n/Na where Na is the length of 
the cubic box. 

The important point is elsewhere: the values of <fi at distinct lattice sites 
are assumed to be independant gaussian random numbers. One can easily 
convince oneself that this implies that the values of the Fourier components if 
and if at distinct sites of the discretized Fourier space are independent gaussian 
random numbers as well. In other words, there is no correlation between modes 
in the initial state. To create a DCC configuration, the quench mechanism has 
not only to be efficient in amplifying the modes (which it is), but it should 
also be able to build correlations between amplified modes. 

2.2 Observables 

Our goal is to determine whether a generic event in the above- described statis- 
tical ensemble looks like an ideal DCC configuration. For this purpose we shall 
compute the event-by-event distribution of the neutral pion fraction f(k) in 

3 In ref. || the variances are given for the lengths of the vectors <f> an d <j>, 
e.g. (0 2 } = E- =1 (^ 2 }^ 2 /4. 



the mode k (see Eq. (||)) on the one hand, and on the other hand we shall mea- 
sure the correlations of this quantity between different modes (see Eq. (|TID). 
Let us briefly review the main ideas underlying the physics we want to describe 
and argue for the relevance of the above mentioned observables. 

In this simplified model, although the size of the box is fixed, we have in 
mind a rapidly expanding system (the quench assumption is a drastic idealiza- 
tion of the effect of expansion). This means that after some time it becomes 
so dilute that the modes decouple and evolve freely. To model this we stop 
the evolution at some freeze-out time tff\. The field configuration at tf is the 
"initial condition" for the subsequent free evolution and determines completely 
the properties of the outgoing waves which describe free propagating pions. 
Indeed, one might consider the classical field configuration at tf as the expec- 
tation value of the corresponding quantum field in the coherent state(we only 
consider the pion sector) 



\a,t f ) = exp I / d 6 k ot(k,t f ) ■ a\k) J |0) , (2) 

with 

3 

a. ■ a' = 2_. a j a ] ) 

3=1 

where aAk) is the creation operator of a free pion with isospin j and mo- 
mentum k, while otj(k,tf) is the eigenvalue of the corresponding annihilation 
operator (tf is a parameter). The state (§) is related to the field configuration 
at freeze-out through (u k = ^k 2 + m\) 

^ tf) = ^MA±^MA. (3) 

\/2uJk 

The subsequent evolution readsf] 

a(k, t>t f ) = a(k, t f ) e~ iu}k ( *-^ . (4) 

In the following, we drop the implicit dependance on the parameter tf, for 
example, oc(k) stands for oc(k,tf) and a(k,t) for a(k,t > tf). Moreover, 
as long as we focus our attention on one particular mode, we shall omit the 



4 A more realistic description would take into account the possible finite time extent of 
the freeze-out period. This would require a model for freeze-out which is out of the scope 
of this paper. We assume that each mode k freezes out at tf which can be thought as the 
end of a common freeze-out period. 

5 The "subsequent evolution" is introduced for the purpose of the argument but is not 
essential. 



indice k, it will be reintroduced when needed. The mean number of quanta 
associated with the classical wave (f|) reads 

fij = (a\a}aj\a) = \aj{i)\ 2 = \ocj\ 2 (5) 

and is time-independant. The fraction / of neutral pions in the mode k is 

/= ~~ • (6) 

Ux + n 2 + U 3 

Let us now discuss the polarisation of the ougoing waves (Q). First, it 
is easy to see that the motion of the vector (p out (t) = ip(t > tf) in isospin 
space is planar. Indeed, the vector I = ip out x <p out is time-independent^. The 
trajectory described by cp out is an ellipse in the plane perpendicular to I. Let 
us call u and L (resp. v and I) the direction and the half-length of the big axe 
(resp. of the small axe) of this ellipse (u 2 = v 2 = 1, u.v = 0, I = ulL u x v). 
One has then 



a(t) = J-(Lu + iIv) e -w-t,mJ ; ( 7 ) 

where rj is some phase factor to be determined from (|J). It follows that 

1x1 I t2 2 i r2 2\ 

n 3 = 2 y L u i +l v i ) ' 
and 

L 2 u 3 2 + l 2 v 3 2 



f 



L 2 + P 



We concentrate on the case of interest, namely the linear polarisation, 
and make explicit the k dependance. A linearly polarised wave is characterised 
by the fact that 1^ = 0, that is: l^ = 0. Then (8^ is the angle between u^ and 
the 7T3-axis) 

a Unear (k,t) = a{k)e' iuJkit -^u n , 

flinear(k) = COS 2 0£ . (9) 

Because both the dynamics (Eq. (P) and the initial ensemble are invariants 
under isospin rotations, there is no privileged direction in isospin space, so 
if such a linearly polarised wave is generically produced, the event-by-event 
distribution of f(k) will be given by the famous l/v7 l aw - This will remain 

6 Ig is the k component of the conserved isorotation generators J d 3 x 4>(x, t) x <j)(x, t) = 
J d 3 kip(k,t) x (fi(k,t). For t > tf the modes are decoupled so each Ij is conserved. 



approximately true for the more realistic case of "flat" elliptic waves (l^ <C Lj~), 
but not otherwise^. The event-by-event distribution of the neutral fraction in 
one mode f(k) is then a useful observable and gives a non trivial information 
about the polarisation of the generic outgoing waves <p^ at . 
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in the same direction 



As emphasized in the introduction, in an ideal DCC all y?2" s oscillate 



So, defining the total neutral pion fraction as 



N3 



pot = £ / 1Q x 

where Nj = j d 3 kfij(k), one has 

fnhc = cos 2 9 , 

where 9 is the angle between u and the 7T3-axis. 

For an ideal DCC, the total neutral pion fraction is distributed accord- 
ing to the 1/y/J law. In a realistic situation, only those modes undergoing 
amplification should be taken into account. Moreover the individual waves 
are neither strictly linearly polarised nor strictly aligned with each others, so 
one would expect some deviation from the 1/y/f law in the distribution of the 
total neutral pion fraction f tot (see e.g. || PI). We will come back to this 
later. In the present work, we compute the following normalized correlation 
function 

c ( %,i) <**>/(?)>. (11) 



U 2 (k))AP(k')) c 

where {A B) c = {A B) — (A) (B) and (...) denotes the average over the ensemble 
at time tf. In the case of linearly polarised individual waves <p°; ut , the /(fc)'s 

measure the directions of polarisation (Eq (|9|)) and C(k,k') is a measure of 
the degree of alignment of these directions. 

7 It is worth noting that, although it is always possible to write f(k) as the squared cosine 
of some angle xti this angle has not, in general, the meaning of a uniform orientation in 
isospin space so f(k) does not, in general, follow the 1/V7 l aw ( see Appendix). In Ref. ||], 
the author has a 1/s/f distribution for f(k = 0), this comes from the fact that the time 
derivative of the field was not included in the definition of the particle number. 

8 The DCC is a zero isospin state: J d 3 xlj: = i. This is very explicitly stated in pF, 
but is also assumed, more or less explicitly, in all the original papers, where the DCC idea 
has been put forward. 



3 Results 

Our box is a 64 3 cubic lattice and the equations of motion (JJ) are implemented 
via a straggered leap-frog scheme |L4| with time-step At = 0.04a. Once an ini- 



tial condition has been chosen, one can follow the time evolution of the <p(k, t)s. 
We reproduce completely the results of Ref. |J with the same bin-average 
procedure in momentum space: in the isospin directions (j = 1, 2, 3) the soft 
modes experience dramatic amplification, the softer the mode the larger the 
amplification, and exhibit coherent oscillations with period 2ir/uJk ~ 27r/m 7r . 
Such behavior does not show up in the a direction. This average behavior^ can 
be qualitatively understood in the mean-field approximation || |15| : for short 
times t < 10a, the curvature of the effective potential (effective mass squared) 
is negative and the softest modes, for which the associated effective frequency 
is imaginary, experience amplification: this is the spinodal instability; for 
times t > 10a, although this phenomenon is no more efficient (the effective 
mass squared is positive), the field is still significantly amplified. This is due 
to the mechanism of parametric resonance triggered by the regular oscillations 
of the mean- field around its assymptotic (vacuum) value |L5|, [16| . Obviously, 



the amplification does not grow forever and at late times (t > 100a) the energy 
is equally dissipated among the modes. In the following we present results for 
two values of the freeze-out time: tf = 10a (21 x 10 3 MC events), correspond- 
ing to the end of the spinodal period, and tf = 56a (10.9 x 10 3 MC events), 
which is the time when the averaged amplification of the soft modes (eq. (|i~2D ) 
is maximal. For simplicity we only sketch results for k = (k = nAk,0,0), 
where Ak = ir/Na ~ 10 MeV (N = 64 is the size of the box). Moreover, the 
linear a-model being an effective theory for energy scales < 100 MeV, we only 
consider the window < n < 15. 

The amplification factor in the mode k at time tf is defined as 

«^=mi- (12 > 

where 

3 

P{k,t f ) =uj k ^2 %(M/) 
i=i 

is the power spectrum in the mode k at freeze-out. The mean numbers fij(k, tf) 
are extracted from the field configuration at time t with the help of Eqs. (S) 



The amplification of a particular mode exhibits large fluctuations around its average 
value over the ensemble. In particular it is not a smooth decreasing function of ||fc|| as 
revealed by an event-by-event analysis. In this respect, the binning procedure of Ref. U 
produces the same result as the ensemble average. The same is observed for the coherent 
behavior of the soft modes. 



and ([|). The fc-dependance of the ensemble average amplification (A(k, tf)) 
is shown in Fig. |l|. It exhibits the expected shape: low momenta are amplified 
with respect to larger ones. Note the semi-quantitative agreement with mean- 
field argument for tf = 10a: from the mechanism of spinodal instability one 
expects a monotonic enhancement of the modes with momenta k < /„-. At tf = 
56a the softest modes have grown further and the window of amplified modes 
has shrunk considerably. 

Event-by-event distributions of the neutral fraction in different modes are 
shown for tf = 10a and tf = 56a in Figs. |3| and [| respectively. All show large 
fluctuations around the mean- value (/) = 1/3. The most amplified modes, 
those which we are interested in, exhibit a very different shape from the less 
amplified ones. A closer analysis reveals that the /-distributions of amplified 
modes are almost the same as in the initial ensemble (see ( |A.5| )), whereas 
those of less amplified modes look very different and all follow the same linear 
behavior. It is interesting to note that this linear law is what is found in a 
thermal ensemble of pions (see ( |A.8[ )), supporting the idea that these modes 
are already thermalised. 

Fig. [5] gives a more precise picture of the above statements. It represents 
the average eccentricity^ of the ellipse (|7|) as a function of the momentum. 
One sees that for tf = 10a the average polarisation is close to its initial value 
for almost the whole studied range of momenta, whereas for tf = 56a only 
the few first modes remain close to their initial value, the average eccentricity 
being constant for k > 50 MeV. Note the correlation between the shapes of the 
average amplification and of the average eccentricity as functions of k. Note 
also that, although the eccentricity of the softest modes decreases as these 
modes are amplified, this is a very tiny effectf] and in fact the shapes of their 
neutral fraction distribution at freeze-out time are mostly determined by the 
corresponding initial distributions. 

So we see that the modes are far from being strictly linearly polarised 
waves, but this causes only a small a deviation from the ideal 1/y/J law. There 
are still important event-by-event fluctuations of the neutral fraction of pions 
with momentum k, so that from the point of view of an experimentalist the 
deviation is not very important. 

Finally, Fig. [1] shows the correlation function (]TT|): at both freeze out 
times the different modes remain completely uncorrelated. In terms of the av- 
erage coherent behavior described at the beginning of this section, this means 
that although the vectors tp^ corresponding to amplified modes do oscillate co- 

10 The eccentricity of an ellipse is defined as the ratio of the lengths of its small to its big 
axes, i.e., in the notations of (0), ek = h/Lk- 

11 An event-by-event analysis shows that the proportion of events in which the zero- mode 
eccentricity is lower than 0.1 is: 13% in the initial state, 16% at tf = 10a and 18% at tf — 
56a. 



herently, their mean directions of oscillation in isospin space (say, the directions 
of their big axes) are completely random: different modes act as independant 
DCCs. The obvious phenomenological consequence is that the signal is washed 
out when summing up the contributions of different momenta as in Eq (|T0|) . 
This is illustrated in Fig. ^ where we show the distributions coming from the 
contribution of five and ten equivalent modes (that is modes having almost 
the same individual neutral fraction distributions). 

4 Discussion 

4.1 Summary 

Starting from a chirally symmetric thermal ensemble for the field <f> = (tt, a), 
we evolved the system with the zero temperature equation of motion of the 
linear a-model. This quenching of initial (thermal) fluctuations is a very effi- 
cient mechanism for producing large and coherent long wavelength oscillations 
of the pion field at intermediate times |J . 

We performed a detailed analysis of the emerging field configuration 
which shows that nothing significant happens concerning its isospin structure. 
The neutral fraction distribution of the amplified, long wavelength modes is es- 
sentially the same as it was in the initial ensemble. Although not exactly l/\/7, 
the distribution is very broad, which is a relevant point for phenomenology. 
However, we found that distinct modes, and in particular the amplified ones, 
have completely independant polarisations in isospin space. This has the im- 
portant phenomenological consequence: the large event-by-event fluctuations 
of the neutral fraction are washed out when the contributions of several modes 
are added, even when one limits one's attention to soft modes only. 

The key point is to realise that the assumption of a completely random- 
ized initial state made in Ref. || implies that, in the initial configuration, the 
field modes are also independant gaussian random numbers. So our result can 
be rephrased as follows: the non-linearity of the dynamics does not build corre- 
lations between modes. The initially postulated chaos is recovered in the final 
state. This contradicts the widely made assumption that the state produced in 
the simplest form of the quench scenario (initially thermalized system) is iden- 
tical to the originally proposed DCC The quenching explains the emergence 
of a strong, coherent pion field, not its hypothetical polarization. 

4.2 Comparison with other works 

Analogous results to that shown in Fig. || were obtained before by several 
authors [|l(J [11], [12| . We briefly review them below. 



Gavin, Gocksch and Pisarski argue in [10] that the system is composed of 
"many small, randomly-oriented domains" resulting in a binomial distribution 
of the neutral fraction of the total number of pions. However, as was pointed 
out later by Rajagopal 0, the system cannot be characterised by a single 
length scale and the picture of different domains of the size of the correlation 
length is not adequate. In the present paper we have shown indeed that the 
typical field configuration emerging from the quench scenario is a superposition 
of misaligned waves in momentum space: different modes and not different 
domains act as independant DCCs^- 



In Ref. ||11| , Rajagopal had in mind a disoriented condensate formed 



by the amplified, long wavelength modes "superposed with short wavelength 
noise" . After bining phase space, he computes the bin-by-bin distribution of 
the neutral fraction of pions with \\k\\ < 300 MeV and observes "an admixture 
of a l/v7 distribution". This momentum space picture is more satisfactory. 
Indeed the observation that modes of different types (amplified vs. thermal- 
ized) contribute incoherently and delay the signal is in line with our present 
result. However, we saw that soft modes do not act together to form a coherent 
misaligned condensate, their polarizations are uncorrelated. 

Finally, Randrup, studying DCC observables in [0, computes the distri- 



bution of the neutral fraction of pions with momenta below different cutoffs. 
From this he argues that "each k contributes pions having an independant 
orientation in isospace" . This is exactly what we have found here. However, 
Randrup's analysis is performed in a model with expansion. This means that 
not all the events entering the analysis have experienced amplification. In this 
sense his work is incomplete and therefore not really conclusive on the matter 
we are interested in here. 

Hence, some of the features of the quench scenario have been noticed in 
earlier studies. However, nobody has carried out a detailed statistical analysis 
of the correlation between mode polarizations. In this respect, we believe, our 
work helps to clarify the situation. 

4.3 Speculations 

As already mentioned, the simplest form of the quench scenario is not sufficient 
to produce a DCC configuration. Inclusion of quantum effects at the mean- 
field level or of expansion will probably not change this. Instead, it would 
be instructive to refine the description by including correlations in the initial 

state. 



12 This does not exclude the possibility a priori that many such domains form in heavy 
ion collisions. However, we argued, in a recent work with A. Krzywicki, that the formation 
of many bubbles undergoing strong amplification is not likely 



10 



A possible way of including initial correlations could be to relax the ther- 
mal approximation for soft modes. Indeed, if the quench mechanism appears 
to be quite natural in the context of high energy heavy ion collisions, it is 
not clear whether the assumption that the initial system is fully thermalized 
is justified. In particular, in such small systems, long wavelength modes may 
not have enough time to thermalize, as it is the case in || (see also Fig.[l]). 
Correlations in the initial state could possibly be amplified by the subsequent 
out of equilibrium evolution. Indeed, from mean-field arguments, the spinodal 
amplification is expected to be quite robust against a large class of initial con- 
ditions^- One has then to have a reliable model for the initial state, that is a 
model for the early stages of the collision. 

Finally, if some correlations are present in a realistic description of our 
initial state, it could be that they do not survive the non-linear dynamics. One 
would be left again with an incoherent superposition of disoriented waves in 
the final state. In such a situation, large fluctuations of the neutral fraction 
of pions in a bin of phase space would only be detectable if one can separate 
individual modes as well as possible. The volume of the bubble should be large 
enough to have a good statistics in each bin, and small enough to have as few 
as possible modes in each bin. 

4.4 Conclusion 

We believe that a complete microscopic understanding of DCC formation in 
heavy ion collisions is not achieved yet. It might also be, however, that the 
original picture of DCC is a too far going idealization. Future investigations 
will tell whether it is the original picture or the quench scenario which has 
to be modified. Theoretical developpements as well as experimental data are 
called for. 



Aknowledgment 

I am very grateful to A. Krzywicki, J. P. Blaizot, and M. Joyce for useful 
conversations and comments and for careful reading of the manuscript. 



13 The main restriction being that initial fluctuations should not be too large in order to 
enter the instability region (negative effective mass squared). 
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A Appendix: neutral fraction distribution in 
a gaussian model 

It is instructive to compute the distribution of the neutral fraction / for the 
following class of statistical ensembles: let us consider the (p^s and the 0g's as 
independant gaussian numbers with zero mean value and dispersions a and b 
respectively (a and b are the same for each mode and for each isospin direction). 
We write the probability (c stands for {j, k}) 

Proba ({(p n }, {0£» = l[ P v (<p c ) P^ c ) dy c d<p c , (A.l) 



1 / x 2 

PJx) = . exp 



^ ' ~ v^W F V 2a* 



1 / x 2 



P ^-7W^{-^)- (A - 2) 

The modes being independant, we focus on one particular k and drop the 
indice. From Eq. ([3D we have for each isospin component 

09 = \ — Rea = A cos 7 , 
V u) 



(p = v 2u Ima = uA sin 7 , (A. 3) 

where 7 is defined through a = y/Ee n , and A = ^2n/uj. The probability 
distribution for the amplitude A and the angle 7 is given by 

P AjJ {A,j) = to AP^A cos 7) P^(u A sin 7) . (A.4) 

The probability distribution of the neutral fraction / = A3 2 j^A 2 + A 2 + A 2 ) 
is 



/+00 / 

dxdydzP A (x) P A (y) P A (z) 5 If - 



x 2 + y 2 + z 2 



where 

f2n 



P A (A)= / d 1 P An (A, 1 ). 
Jo 

After some calculations one gets 

P f (f) = \(Fn(f) + F-n(f)), (A.5) 

where 

w) = (»-d-/))( ii _ Si (1 + _ 1 2/) N3/2 

12 



and 

b 2 + ujg? 

b 2 — uoa 2 
It is instructive to look at some specific examples of (|A.2| ). First, we can 



fix one of the two vectors <p and cp to be zero. In this case a is proportional 
to a randomly oriented vector in isospin space and one recovers the l/v7 law. 
Indeed, let us choose b = (Q = —1), that is ip = 0. Then ex. ~ cp and 

W) = ^J. (A.6) 

The second interesting example is the case where b = uoa (Q = +oo). 
In terms of a probability distribution for the complex numbers aj = Xj + iy$^\ 

P a (oi) d 2 a = P x , y (x, y) dxdy , 



one has, 


from 


(|4.2) 


and 


(A.3), 












P«{a) = 


where 


















a 2 


= ab = 


--(n) 



Tea 2 



2 / 2 

\a\ /a 



d aP n {a) \a\ 



(A.7) 



(|a| 2 = n is the mean number of quanta in the coherent state \a) and (n) in 

the mean number of quanta in the statistical ensemble). 

For the class of ensembles ( |A.7| ), the neutral fraction distribution reads 

P f {f) = 2(1 - /) . (A.8) 

A particular case of ( |A.7|) is that of a thermal collection of quanta with fre- 
quencies uj. 

Finally, the initial ensemble used in the quench scenario belongs to the 
class ( |A.1| ). The 0j(x)'s and the 0j(x)'s are independant gaussian numbers 
of variances (<fr 2 ) and (<fi 2 ) respectively. This implies that the (pj(k)'s and 

the 0j(/c)'s are independant gaussian numbers of variances M {(ft 2 ) and M {(f> 2 ) 
respectively [M is a normalisation factor) 

a 2 = Ar($)=Ar^,b 2 = Ar(<p 2 )=Ar v ^, 

and 

A-u 2 
Fig. shows the neutral fraction distribution in the initial ensemble for the 
two extremal values of the studied window < k < 150 MeV0. 

14 Such a probability distribution for the a's is called the P-representation of the den- 
sity operator |17| : a very large class of quantum density operators can be written p = 
Jd 2 aP a (a)\a) (a\. 

15 Note that although our initial state is a thermal collection of quanta, these are not zero 
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temperature pions (this means that their dispersion relations are not those of free pions). 
That is why the initial neutral fraction distribution is not the same as (A. 
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Figure 1: The average amplification (A(k,tf)) as a function of k (MeV) for tf = 
10a (left) and tf = 56a (right). The error bars represent statistical errors and the 
lines are just guides for the eyes. 





Figure 2: Neutral fraction distribution in the initial ensemble for n = (left) and 
n = 15 (right). 
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Figure 3: Neutral fraction distributions at time tt = 10a for n = (left) 
and n = 15 (right) (k = nAk with AA: ~ 10 MeV). All the histograms for n < 14 
are compatible with the corresponding initial distributions while n = 15 already 
exhibit a linear shape to be compared with the 2(1 — /) neutral fraction distribution 



in a thermal ensemble (Eq. (A. 8)). 
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Figure 4: Neutral fraction distributions at time tt = 56a for n = (left) 
and n = 4 (right) (k = nAk with Ak ~ 10 MeV). The histograms for n < 3 
follow the corresponding initial distributions. All the modes n > 4 exhibit the same 
linear shape. 



16 



ft, 



s 



3 
2 

1 



0.2 0.4 0.6 0.8 



./' 



0.2 0.4 0.6 0.8 



./' 



Figure 5: Neutral fraction distributions where we took into account the contribu- 
tions of modes with n = 0, ...,4 (left) and n = 0, ..., 10 (right) at tf = 10a. At this 
time all these modes have essentially the same individual distribution (cf. Fig. |3|). 
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Figure 6: The average excentricity (e(k,tf)) as a function of k (MeV) in the initial 
ensemble (square) and for tf = 10a (diamond) and tf = 56a (triangle). The lines 
are just guides for the eyes. 
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Figure 7: The correlation function C(0, k) (Eq. (|TTJ) ) vs. k (MeV) for tf = 5a (left) 
and tf = 56a (right). We computed also C(ko,k) as a function of k for different 
values of ko, all look the same as above. 



